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Abstract 

We investigate the phenomenon of the decay of a supercurrent due to the 
zero-temperature quantum tunneling of vortices from the edge in a thin su- 
perconducting film in the absence of an external magnetic field. An explicit 
formula is derived for the tunneling rate of vortices, which are subject to the 
Magnus force induced by the supercurrent, through the Coulomb-like poten- 
tial barrier binding them to the film's edge. Our approach ensues from the 
non-relativistic version of a Schwinger-type calculation for the decay of the 2D 
vacuum previously employed for describing vortex-antivortex pair-nucleation 
in the bulk of the sample. In the dissipation-dominated limit, our explicit 
edge-tunneling formula yields numerical estimates which are compared with 
those obtained for bulk-nucleation to show that both mechanisms are possible 
for the decay of a supercurrent. 

I. INTRODUCTION 



Recently there has been much renewed interest in the physics of vortices in 
high-temperature superconductors in the presence of applied magnetic fields 
p). Phenomena involving vortex dynamics and their hindrance by means of 
added pinning centres pose well-defined physics problems of great relevance 
for the technological applications of the new high-T c materials. 

Some of the lesser understood issues are related to the problem of the 
residual resistence due to the thermal or quantum motion of the vortices in 
the presence of an externally applied supercurrent, but in the absence of the 
applied magnetic field. Even in the absence of an external field, vortices can be 
dragged into the bulk under the action of the Lorentz-like Magnus force which 
they experience sitting at the edge of the sample. This phenomenon, leading 
to an electrical resistence in the presence of dissipation, has been investigated 
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recently by relatively few authors |2],||] considering its importance from both 
the conceptual and the practical points of view. In addition, in a recent 
article we have pointed out how such a residual resistence may arise, in 
the absence of an applied field, also due to the spontaneous homogeneous 
nucleation of vortex-antivortex pairs in the bulk of the sample. These are 
created as fluctuations of the electromagnetic-like Magnus field acting on the 
quantised vortices and antivortices thought of as electron-positron-like pairs. 
By means of this analogy with quantum electrodynamics (QED), we have set 
up a powerful "relativistic" quantum field theory (QFT) approach to study 
vortex nucleation in the two-dimensional (2-D) geometry. This study has been 
conducted in the presence of quantum dissipation and for the cases of either 
an harmonic local pinning potential jjj or, more recently ||, for a periodic- 
lattice distribution of pinning centers. Our analysis has produced an explicit 
analytic result, and an estimate of the magnitude of the vortex pair- nucleation 
rate has shown that the effect may become experimentally accessible at low 
temperatures and high enough current densities. 

In this article we investigate, by means of our QFT approach, the origi- 
nal issue of tunneling of quantised vortices from the sample's boundary as an 
alternative mechanism for the decay of an applied supercurrent in a thin-film 
geometry. The problem has already been introduced by Ao and Thouless [|| 
and by Stephen |J, with and without the inclusion of quantum dissipation, 
and discussions that give full weight to the inertial mass of the moving vor- 
tex can be found in the literature alongside treatments [^J6| in which the 
inertial mass is taken to be negligible. The physics of the problem has been 
captured mainly through semi-classical evaluations of the tunneling rate and 
discussion has centered on the effects of dissipation and pinning in opposing 
the stabilising influence of the Magnus force on the classical orbits of a vortex 
H]||. Beside the issue of the residual superconductor's resistence, the quan- 
tum tunneling of the magnetic flux lines into the bulk of the material has 
been investigated to understand theoretically ||[7j the observed || flux-creep 
phenomenon in the presence of a field. 

In this work we will analyze further the problem of vortex tunneling from 
the point of view of the instability of the "vacuum" represented by a thin su- 
perconducting film in the presence of an externally-driven supercurrent. Our 
aim is in fact that of obtaining an explicit formula (complete of prefactors) for 
the tunneling rate. From the theoretical point of view, our approach differs 
from other semi-classical evaluations of the path integral in that use is made of 
the Schwinger formalism for pair creation in QFT. Since this fully-relativistic 
formulation provides for a description of vortex-antivortex pair nucleation, 
we will be concerned here with the completely non-relativistic limit of this 
approach, formally corresponding to the case of a negligible vortex inertial 
mass and to contributions to the path integral relating to quantum particles 
moving "forward" in the time coordinate. This non-relativistic version of the 
Schwinger formalism indeed appears as very convenient for carrying out ex- 
plicit and straightforward calculations, especially in the dissipation-dominated 



2 



regime. We in fact work out the tunneling rate in the case when dissipation 
dominates over inertia, and for the case of a potential barrier made up of the 
Coulomb-like attraction to the edge and the electric-like potential extracting 
the vortex into the bulk. Our central result shows that the tunneling rate T 
has a strong exponential dependence on the number current density J, much 
as in the case of vortex nucleation in the bulk of the sample Q . The numerical 
value of the tunneling rate per unit length, as obtained from our approach 
with material parameters typical of the high-temperature cuprate supercon- 
ductors, compares favourably with analogous results obtained for the bulk 
nucleation rate. This qualifies the novel vortex nucleation mechanism by us 
proposed in recent articles [Qj5| as almost as likely as vortex tunneling from 
the sample's edge for the decay of a supercurrent. In our estimates, edge 
tunneling appears to be more favourable than bulk nucleation for equal given 
extension of edge length and surface area. 

It is useful to describe the problem at hand by means of the classical 
equation of motion for a single vortex moving at relatively low velocity in a 
supercurrent having density J 

mq = — VJ7(q) + eE — eqxB — r/q (1.1) 

Here m is the (negligible) inertial mass of the vortex carrying topological 
charge e = ±2n and treated as a single, point-like particle of 2-D coordinate 
q(i). Also, U(q) is the phenomenological potential acting on the vortex, in 
the present case the "Coulomb-like" interaction binding it to the sample's 
edge. The supercurrent J gives rise to an electric- like field E = xJ (a no- 

(3) 

tation implying E • J = 0) superposed to a magnetic-like field B = zdp s 

for a thin film orthogonal to the vector z and having thickness d with a su- 

(3) 

perfluid component characterised by a 3-D number density p s . Finally, 77 
is a phenomenological friction coefficient taking dissipation into account. A 
derivation of this electromagnetic analogy for the Magnus force was given by 
us in Ref. pj (but see also pHIII]). Further contributions to B arising from 



other quantum effects are possible and have been recently debated [12|, how- 
ever in the dissipation-dominated regime the actual explicit dependence of B 
on the material's parameters will be irrelevant for our study. The quantum- 
mechanical counterpart of Eq. ( |1.1| ) is constructed through the Feynman 
path-integral transposition in which the dissipation is treated quantistically 



through the formulation due to Caldeira and Leggett [13|. This approach 
views quantum dissipation as described by the linear coupling of the vortex 
coordinate to the coordinates of a bath of harmonic oscillators of prescribed 
dynamics. By means of this description of quantum dissipation, we formulate 
our own approach to the dissipative tunneling of vortices subject to the fields 
E and B (the magnetic-like field being treated in an effective way). 

The organisation of our article is as follows. In Section II we sketch some 
basic facts about our Schwinger method for the relativistic vortex dynamics 
and show how the completely non-relativistic limit yields a convenient formu- 
lation of the vortex tunneling problem. This is expanded in Section III, where 
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we treat in detail the situation in which dissipation and a Coulomb interac- 
tion are present and dominate over vortex inertia. In Section IV we estimate 
the tunneling rate for a range of reasonable values of the material parameters 
and compare the results with previous estimates for the bulk nucleation rate. 
This Section contains also our conclusions. We work in the units system for 
which H = 1. 



II. VORTEX QUANTUM DYNAMICS: A UNIFIED APPROACH TO VACUUM 
DECAY THROUGH NUCLEATION AND TUNNELING 

The calculation that follows is the entirely non-relativistic version of our 
published QFT formulation GJ for the decay, via vortex-antivortex pair nucle- 
ation, of the "vacuum" represented by a supercurrent having number density 
J flowing in a superconducting thin film. Since the would-be particles should 
experience a force entirely analogous to the Lorentz force from a uniform elec- 
tromagnetic field, we expect this vacuum to be unstable and decay via vortex 
pair creation. We evaluate the probability amplitude for the vacuum decay 
in time T, 

Z = (0\e- iAr \0)=e- iTWo (2.1) 

where Wq = £(vac) — i^ gives the energy of this vacuum, £(vac), and its decay 
rate, T. With a suitable normalization factor, the probability amplitude is 
given by a functional integral over field configurations. In the presence of a 
gauge field A^(r,t) and external potential V(r), we have 

Z = M J Vcj)exp |-*y d 2 r dt 0* \ —Dq + ^D 2 - £ 2 - V(r)J 

= exp | -Tr In ^— — D 2 - Z>| + £ 2 + F(r)J j (2.2) 

which can be evaluated, in the Euclidean metric, by means of the identity 
-D% + £l + V{v) 



Tr In- 



A 2 



lim/ ^Tr{e-(-^o 2 +V>_ e -AVj (23) 



is 



where (with x a =it) D 2 E = D 2 + ^{D\ + D%) and where = <9 M - iA^ 
the usual covariant derivative. Notice that 7 = m/£o, the ratio between the 
inertial mass m and the activation energy £q, plays the role of the inverse 
square of the velocity of light. Thus the non-relativistic limit is implicitly 
taken in the dissipation-dominated case, where m^O. 

The evaluation of the trace is straightforward by means of the Feynman 
path-integral. We have, with a suitable normalization factor M(t) 

Tr{e- ( - D z +£ Z +v)T \=Af(T) [ dq [ Pg(s)e _ X ' dsCB (2.4) 

L J J Jg(0)=q(r)=g o 



4 



where the Euclidean version of the relativistic Lagrangian reads 



Here q is taken to be a (d + l)-dimensional coordinate describing the closed 
trajectories in the Euclidean space-time as a function of the Schwinger proper 
time s. We denote q = (qt(s) , q(s)) in terms of its time- and space-like 
components, respectively, the dot representing the derivative <? = ■ The 
part of the trajectory moving backward in time qt is therefore interpreted 
as describing the antivortex. The expression for the vacuum decay rate is 
therefore, quite generally 



L=l m [°°— f (iqoAA(r) / Vq(s)e- ^ dsCs 

2 Jo T J Ja(G)=aM=an 



q(0)=q( T )=q 



where the integral over qo, the initial point in space of the closed paths, runs 
over the region of particle nucleation. In the case of vortex-antivortex pair 
creation in the bulk of the thin film, this region is the surface of area I? of 
the film. 

In the problem considered in the present paper, vortices are already 
present in a narrow strip of width a~£ (with £ the coherence length of the 
superconducting order parameter) along the edge of the sample where V x J 
is different from zero. Thus the "vacuum" in this case corresponds to a uni- 
form distribution of vortices along the edge, whilst none of them is present in 
the bulk in the absence of an external magnetic field. The "vacuum decay" 
thus corresponds to the possibility that some of the vortices are dragged away 
from the edge and enter the bulk. We can still take Eq. fl2.6|) as the starting 
point for determining the decay rate, provided we factor out the contributions 
coming from all paths describing particles propagating backwards in time qt 
and replace them with a suitable normalization factor. Also, since vortices are 
already present, a chemical potential /i* will be introduced to cancel out their 
nucleation energy. We will take the non-relativistic limit explicitely, since now 
the advantage of the relativistic formulation (which in our case provided for 
a simultaneous description of both particles and antiparticles) is no longer 
useful. 

We find that this way of approaching the tunneling problem, based on Eq. 
( |2.6j ), has in our opinion a number of advantages over the standard instanton- 
type calculation (7,14|. The instanton (or WKB for point-like particles) cal- 



culation, although well established, nevertheless calls for the determination 
of the "bounce" solution of the path integral's saddle point equation and the 
functional integration of the fluctuations around it, a task often too hard 
to carry out explicitely. The present formulation based on an entirely non- 
relativistic Schwinger approach leads to a promising viable alternative, as we 
now show. The integration over the time paths qt(s) in Eq. ( p.q 1 ) is carried 
out by means of a sadde-point approximation which fixes the correct non- 
relativistic form of the Lagrangian. Our uniform-field situation corresponds 



to A^(q) = \F iH ,q v , with the (2+l)-dimensional field tensor given by, after 
the analytic continuation to Euclidean time 




B iE x 

Ffiv = | —B iE y | (2.7) 
-iE y 

In the present treatment, the vortex chemical potential fj,* is also added to 
the scalar potential, Aq — > Aq + fj,*, so that with a suitable choice the forward- 
moving vortices can be selected from the backward-moving antivortices. In 
our previous work Q we have shown that the main role of the magnetic- 
like field in the nucleation of vortex-antivortex pairs is to renormalise the 
friction coefficient (mimicked by the bath of harmonic oscillators included in 
the potential V(q)), 

r}^m = {r, 2 + B 2 )h (2.8) 

(denoted simply by r\ in what follows), as well as the nucleation energy 

£ 2 ^£ 2 R = £ 2 + A£(B) 2 (2.9) 

where A£(B) is a B- and 77-dependent renormalization. However, the friction 
coefficient 77 is a rather uncertain normal-metal parameter and, moreover, we 
assume that there is an infinite reservoir of vortices at the border of the film, 
their /i* cancelling the value of the effective nucleation energy as well as all its 
renormalizations. Thus, we may safely neglect the effects of the magnetic-like 
component of the Magnus force and write the relativistic action in the form 

S E (B = 0) = £ ds ||q 2 + ^q 2 - q t (B ■ q + M *) + £ 2 m + F(q)| (2.10) 

= £ ds^q 2 + \ (q t - 2(E • q + //)) 2 - (E • q + //) 2 + £ 2 R + F(q)} 

The saddle point approximation fixes the function qt = dqt/ds so as to have 
an extremum for the action: namely, we have the saddle-point condition 

^|P = ife-2(E.q + ^))=0 (2.11) 

from which we obtain 

^ =2(E-q + /x*) ^2^ (2.12) 
as 

where we take /i* = ±£qr to be the dominant energy scale. Here, the positive 
sign applies to particles moving forward in the time qt and the negative one 
to backward-moving antiparticles. The resulting saddle-point action for those 
trajectories moving forward up to the standard Euclidean time T = 2£q R t 
can therefore be taken as (indicating for short with t the time component qt) 
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5+(0)« / dt{-ml-£) -E-d + U(q)>=SNR (2.13) 



where the residual contribution in terms of the nucleation energy is cancelled 
by our choice for the chemical potential and where a term (E • q) 2 /2£q R 
has been dropped. Also, we have denoted by U(q) = V '(q) /2£q R the non- 
relativistic potential energy. We now carry out the functional integral over 
the trajectories moving forward in time, replacing the contribution of the 
integral over the antivortex trajectories with some factor <]?. This factor will 
also include the contribution arising from the fluctuations around the saddle 
point. The resulting expression for the tunneling rate through the barrier 
represented by the overall potential U(q) — E • q becomes, from Eq. 

-dr.., _ r _ - Snr 



-«/m / —$M NR (T) / 2?q(f)e- 4Nfl (2.14) 

2 JO J- Jq(T)=q(0) 

This expression contains only closed paths forward-moving in time t and 
weighted by the non-relativistic action. The standard normalization factor 
Mnr(T) corresponds to the path integral for the non-relativistic free particle 
in D = 2 dimensions, in the absence of dissipation 

N NR (T) [ p q( t )e -/o T ^W = ™ (215) 

7q(0)=q(T) 2?rT 

To fix the overall factor we proceed in the following manner. The expression 
for the non-relativistic vacuum decay amplitude is as follows 

(0\e- HNRT \0) =AT NR (T) J Vc[fT SNR 

= exp - (£{vac) - i^j T = e~ £ ^ T (cos ^ + i sin ^ (2.16) 

Since we are dealing with a highly-stable vacuum, we take the lowest order 
term in the expansion in powers of TT for T<^.£(v ac) to obtain 

Vqe-SNR = i_ e -£(vac)T + . . . ( 2 J7 ) 

which determines, inserted in Eq. ( |2.14j ), the following expression for <I> 

$ = £{vac) = ( f°° dTe- £i - vac ^ T \ (2.18) 



o 

In deriving this expression we have taken into account the fact that the cal- 
culation for the rate T is based on a saddle-point approximation where $ is 
a function of the material's parameters only. The exponential e - £ ( vac ) T rep- 
resents the vacuum probability amplitude in the absence of the electric field 
instability, thus we finally write 

r _ Im Jo 00 ^M NR {T) J Vqe- s ™ 
2 ffdTJVNRWfVqe-!®* 



(2.19) 
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where the action Sj^ R is the non-relativistic action in the absence of the elec- 
tric field. We have in this way reached a formulation for the tunneling rate 
which leads to an expression in agreement with the standard WKB method. 
This can be checked by examples insofar as the exponential term (3 is con- 
cerned, in the characteristic expression for the tunneling rate T = ae~^ ' . 



III. EVALUATION OF THE TUNNELING RATE THROUGH THE EDGE 

POTENTIAL. 

In this Section we will apply the above general formalism to the situation 
at hand, where the vortices tunnel from the edge strip under the combined 
effect of the "electric potential" — Eq and the 2-D "Coulomb electrostatic po- 
tential" . The latter is due to the attraction of the vortex by the edge, which is 
equivalent to the attraction by a virtual antivortex, like in the familiar virtual- 
charge method in standard electrostatics. This Coulomb potential takes one 
of the equivalent forms Uc{y) = K ln(l + y/a) or Uc{y) = \K\a (1 + (y/a) 2 ), 
where y is the distance from the edge, both acceptable extrapolations of the 
known large-distance behavior. The coupling constant K depends on the 
superfluid density pf^ and carrier's mass rriQ through the relationship p| 
K = d/rriQ. However, its precise value in real materials is unknown 

and this parameter will be treated, like many others in this calculation, phe- 
nomenologically. 

When the Caldeira-Leggett quantum dissipation is taken into account 
[13,|], we end up with an overall non-relativistic potential 



U(d) = ^-V(q) = U c (y) + U D (d) 
^(q) = 53 | \ m ^k + \ m k^l (xfc + 



V 

2m k ujl I 



with the oscillators' masses m/% and frequencies u) k constrained in such a way 
that the classical equation of motion for the resulting action reproduces the 



(3.1) 



E ~ w*)= J(w) = rpj (3.2) 

' ^ rri k uj k 



form (1.1). This requires the c k to satisfy the constraint (for the so-called 
ohmic case, which we consider) 

9 

IT ■ 
~ k 

rj being the phenomenological friction coefficient of Eq. ( |1.1| ) (renormalised 
by B). Notice that this coefficient is unaffected by the non-relativistic limit. 
After a Fourier transformation in which the x& modes can be integrated out 
Q , we are lead to an effective non-relativistic action which in the dissipation- 
dominated (m — > 0) regime reads, with q = (x, y) 

S NR = 27T// V nq n -ql-ETy + K f dt In ( 1 + ^) (3.3) 

Jo \ a J 
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Here q(i) = q + J2 n ^0 ( In e * Wn *j where to n = 2irn/T, and furthermore y = 
Jq dty(t) /T stands for the n = mode in the y-direction. In order to render 
the problem tractable, at this point, we introduce a further approximation for 
the interaction term of the action in Eq. (p.3|) , by writing 



T dtln (l + -) «Tln ( 1 + - 



which is justified in the limit of large exit times T (see Appendix). We are 
therefore in a position to evaluate the numerator of the formula, Eq. (2.19), 
for the tunneling rate, which we write schematically as T/2 = M I'D. 
We define, for the numerator 

poo (VT f roo fVT 

M=Im J —Mnr(T) J V^e~ SNR = Im J —M NR {T)I x I y 

where the path integral over q(£) factorises into I x = L Y\'^L 1 (l/2r]n), repre- 
senting the contribution of the free dissipative motion along the x-axis of the 
edge having length L, and in 



ly 



J dyI G (y,T)e ET y- KTl ^ 1+ y^ 



the factor containing the dynamical effects, which we now evaluate. This we 
do most conveniently by integrating out first the real and imaginary parts, 
Re y n = ip n and Im y n = £ n , of the modes y n for n > 0. Since the path 
integral is restricted to loops which have their origin at y(0) = y{T)=yo = 0, 
the constraint yo = y + Sn^oVVi = V + 2E(i>o^ii = has to be imposed, 
corresponding to the vacuum wave-function | ^0(2/0) 1 (2/0) - This leads to 
the constrained Gaussian integrals 

„ 00 / 00 \ -. 

I G (y, T) = JH ly + 2 £ VU e^E^ = _ e ~W x 

n=l \ n=l / " 

x Yl ( g — ) / II d ^ n ex P { ~ 27Tr l ^n(nS nm + l)ip m - 2nr]y ^ Vv, 

n =l \ Z V n / J n=2 I n,m=2 n=2 



where the 5- function constraint has been taken care of through, e.g., the ip\- 
integration. Introducing the matrix (n, m = 2, 3, 4. . .) 

M nm = 2-KT](n5 nm + 1) 

the integrals can be evaluated, formally, to give 



m=2 "^nm. 



Notice that the resulting expression involves infinite products and sums that 
would lead to divergencies. However, since the dissipation is believed to 
be suppressed above a characteristic material-dependent frequency u ci these 



products and sums are to be cut at a characteristic integer n* = [u c T/2ir]. 
Writing M = 2ttt]M (1 + M^V), with M 0nm = nS nm and V nm = 1, we get 

n* 

detM= JJ(2vrr/n) exp Tr ln(l + M^V) (3.10) 



n=2 



The trace can be evaluated now by formally expanding ln(l + M Q V), to get 

n* , 

det(l + M^V) = 1 + ~=Q ( 3 - n ) 

n=2 U 

In a similar way we evaluate = ^(^ nm - q^)> yielding 

«* 1 -. 

S^-^-q' (3 ' 12) 

n,m=2 1 ^ 

from which we obtain 



2\QJ ^\2 V U 

The numerator of our formula for the tunneling rate T is therefore the double 
integral 



^=^K?) 1/2 rf«n(^)7;*- ss ~ i " (,+s/ " ) <-> 



We remark that in the above formula the following formal expression appears: 
N=Mnr{T) ]ln>o (^n) • interpret this expression as an m —> limit: 



-►— exp-TAf (3.15) 

27T 

where AS = — (1 + ln(ma; c /r/)) is a further renormalization of the activation 
energy which is also cancelled by tuning the vortex chemical potential fi*, so 
that, effectively, N—y-^. 

The integral is now evaluated by means of the steepest descent approximation 
in both y and T; the saddle-point equations are, with 

S=^y 2 -ETy + KTln(l + ^) (3.16) 

as the reduced action, 

dS v 

— = -Ey + KHl + ^) = 

^ = ^y-ET + KT^-_ = (3.17) 
dy Q a + y 
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the first of which yields y and the second T, the mean exit distance and 
time, respectively (since n* is a large cutoff integer, Q is kept fixed in the 
variation of the effective action). This leads to the expression Sq = -K7]y 2 /2Q 
for the saddle-point effective action. One has, furthermore, to integrate over 
the Gaussian fluctuations around y and T. One can easily verify that there 
is a negative mode, thus giving an imaginary contribution, coming from the 
integration over the fluctuations in T. We in fact expand 

S = S + a5y 2 + QSyST = S + a{5y + — 5T) 2 - ^-5T 2 (3.18) 

2a 4a 

with a = — tf^gy-i ) aim (3 = —E + This leads to the following 



evaluation of the numerator of Eq. ( p. 19 ) 



dr., ,„,/■_ - nr e~ So 



roo fTT p 

J\f = Im —N NR {T) J VqeT SNR 



o T JV " V 'J " ' V 2QT(E-K/(a + y)) 



S = ^f (3.19) 

where Q = ]n[uj c T /2ir] for large enough T. 

We now come to the evaluation of the normalization denominator in Eq. 
(2.19). We remark that, the numerator being already of the form ae~@, this 



denominator contributes only to the specification of the prefactor a in the 
expression for the tunneling rate T. To begin with, we carry out the Gaussian 
integrals over £ n and ip n . Here we use the alternative form of the Coulomb 
interaction, Uq = \Khx (l + (y(t)/a) 2 ), and we expand the logarithm around 
the minimum (now at y = 0) to lowest order (thus keeping all the integrals 
Gaussian). We define 

pOO p (q\ poo 

V= J dTM NR {T) J V^mr = J dTN NR (T)I x I^ (3.20) 
where I x is the same as for the numerator above, while 

OO „ CO 

W d£ n dip n / dy5(y + 2^2 Vv, 



r(o) 
y 



oo 

'>n) X 

n=l '' n=l 

{oo fcCT 00 1 

-27TT, J2 <tl + €) -^EKn + + f\ \ (3-21) 

n=l n=l ) 

Finally, repeating many of the steps of the calculation done for the Gaussian 
integrals of the numerator J\f, we end up with 



^2^2 poo ™ n 2 n 2 i ; 

V = —J—L / du (1 + 2uQ(u)y 1/2 TT (1 + u/n)' 1 =^L ^a 2 r,) (3.22) 

K Jo ^ K K 

,2, 



where u = KT/2ira r], Q{u) = J2n=i^/i n + u ) an d <P * s a dimensionless 
function of order unity. The evaluation of the tunneling rate is now formally 
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completed, and an explicit analytic formula ensues from our (albeit approx- 
imate) treatment leading, from Eq.s ( |3.19f ) and (|3.22|) , to the tunneling rate 
per vortex: 

T K e~ s ° 

-=N/V= -^==_ _ _ - — (3.23) 

This expression is readily evaluated, once the parameters a,r],u c ,K and J are 
fixed, by evaluating numerically the solutions of Eq. ( |3.17| ) and the integral 
in Eq. (|3.22|) by means of n* = [uj c T /2tt] = [(uj c a 2 r]/K)u] ([] denoting the 
integer part). To obtain the tunneling rate per unit length, we must compute 
the density of vortices along the edge. Since p v = x V<9 is the bulk 
density in terms of the condensate phase 9, and since the number current is 

(3) 7 

J = ^{S/e - A), we obtain 

/ dyp v = -^L- [ dyd y J x = J/K (3.24) 
The final expression for R, the tunneling rate per unit length, is thus 

R = ^T (3.25) 

We conclude this Section by noticing that, since E = 2ttJ in terms of the 
supercurrent density, and since the saddle-point equation ( 3.171 ) has iterative 



solution y = ln(l+y/a) = ^{ln(l+K/Ea) + - ■ •}, we can cast the dominant 
exponential factor determining T in the form 

r]{KHl + K/2irJa) + ---} 2 
S ° = 8^PQ (3 ' 26) 

This yields almost the same leading dependence (logarithmic corrections 
apart) T~e~( J °/ J,> on the external current that was also obtained for the ho- 
mogeneous bulk-nucleation phenomenon |3j (in terms of the tunneling length 
£ T r^K/E we could also rewrite r~e -r ^r). Indeed, the effective Coulomb en- 
ergy K ln(l + K/2tt J a) can be interpreted as playing the role of the activation 
energy £q renormalised by the vortex- antivortex interactions. 



IV. DISCUSSION AND ESTIMATE OF THE TUNNELING RATE. 

We have presented a new treatment for the quantum tunneling of vor- 
tices from the boundary of a thin superconducting film. The formulation 
makes use of the idea that tunneling can be viewed as a pair-creation process 
in which only forward-moving time trajectories contribute to the Schwinger 
path integral. Backward-moving antivortices represent contributions to the 
path integral that are factorised out, and we have made use of the fact that an 
infinite reservoir of vortices is present at the edge. This leads to a mathemat- 
ically convenient formulation of the evaluation of the tunneling rate, which 
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we have worked out in detail by assuming that the vortices experience an at- 
tractive 2D Coulomb-like potential confining them to the edge. By means of 
some approximate treatment of the logarithmic Coulomb interaction, and of 
the saddle-point evaluation of the path integral, we have reduced all integra- 
tions to Gaussian integrals which afford a closed analytical expression for the 
tunneling rate. This can then be estimated by chosing material parameters 
suitable for typical cuprate superconducting systems, e.g. YBCO films, that 
afford some of the highest critical current densities ( J c ~10 7 Acm~ 2 at 77 K 
[15]). To fix the indicative value of the friction coefficient rj, we make use of 
the the Bardeen-Stephen formula [16] linking the upper critical field B C 2 to 
the normal metal resistivity p n : r] = <&qB C 2/ p n ($0 being the flux quantum). 
This yields r/~10 -2 A~ 2 for YBCO. A rather uncertain parameter is the cutoff 
frequency u c , although the dependence on it of the tunneling rate is rather 
weak. We therefore take indicatively w c «80 K, just under the expected value 
of the Kosterlitz-Thouless transition temperature in these films Q . As for the 
edge thickness, we take a~£~10A of the order of magnitude of the coherence 
length. Finally, there is the value of the Coulomb interaction strength, K. 
Since we must have K < uj c , we present in Fig. 1 the dependence of the tun- 
neling rate on the current density for some indicative values of K. The effect 
of temperature is taken into account phenomenologically Q, by introducing 
thermal current fluctuations in the exponential of the formula for the rate, 
J 2 -» J 2 + AJ(T) 2 . The fluctuation term AJ(T) is fixed by requiring that 
for J = an Arrhenius form T cx e~ Umax ^ T is obtained for the rate, with the 
barrier's height U max = K ln(K/2irJa) — K + 2tt Ja as the activation energy. 
Notice that the barrier is infinite for J^O, so that there is no zero-current 
tunneling of vortices even for non-zero temperature. It can be seen from Fig. 
1 that the effect of thermal fluctuations, as estimated from the above inter- 
polative assumption, is indicatively to raise the tunneling rate by some orders 
of magnitude. 

The results are summarised in Fig. 1, where it is shown that tunneling of 
vortices from the film's boundary can become a highly likely phenomenon for 
sufficiently large current densities J and relatively low Coulomb couplings K. 
The process is aided by the presence of thermal current fluctuations. For even 
larger currents than those considered here, the process could be described by 
means of a classical treatment based on the time-dependent Ginzburg-Landau 
equation p7[j. It is interesting to compare the results obtained in Fig. 1 with 
the estimates carried out for the bulk nucleation process Q. There, the 
maximum nucleation rates observed (also taking temperature into account) 
were in the region of ~ 10 13 ^m _2 s _1 for a current J = 10 7 A cm -2 and 
assuming that pinning centers act as nucleation seeds. We therefore conclude 
that edge tunneling appears to be a more likely supercurrent decay mechanism 
for normal superconducting sample geometries. Yet, the bulk- nucleation of 
vort ex- ant i vortex pairs is also a competitive parallel mechanism that could 
be singled-out by means of an appropriate choice of sample geometry. 

Although the measurements of the residual resistence in a superconductor 



13 



may not allow to distinguish between these two competing mechanisms, we 
stress that from the point of view of vortex-tunneling microscopy the two pro- 
cesses remain entirely separate observable phenomena. Bulk pair-nucleation 
remains, in particular, a completely open challenge for observation of a phe- 
nomenon that even in the context of standard QED has remained, to our 
knowledge, so far elusive for static fields. 
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APPENDIX 



In this Appendix, we justify the use of the approximation, Eq. (3.4), used 
in our calculation for handling the "Coulomb" interaction in the path integral 
of the numerator of our formula for the nucleation rate. For large y, and with 
<jj n = 2Trn/T: 



o \ a 



Tin l + £ +YJ 



1=1 



1 



dt y ^ Wnt 



At the saddle point, y n = —y/2nQ (as can be easily verified by varying Eq. 
(|3.3|)) and therefore we have, for example 



[ T dt (E 



Vn c iu n t 

y 



4Q< 



E 

n=l 



n- 



Since Q = ln[^-] and T~^-, we get the estimate 

; T 



o V a 



-2' 



Thus the approximation of taking the "Coulomb" potential of the form 
K ln(y/a) is justified for very low currents (E = 2ttJ— >0). We still keep this 
form, when the current is not so low, as we believe it captures the essential 
point of the dynamical picture. 



(A.l) 



(A.2) 



(A.3) 
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FIGURE CAPTIONS 

Figure 1. Plot of the tunneling rate per unit length, R, as function of the 
supercurrent density, J, for the values K = 30 K and K = 50 K. Dashed 
lines, T = K; full lines, T = 2.5 K. 



15 



